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A direct fit of the isoscalar spin-orbit and both isoscalar and isovector tensor coupling 
constants to the / 5 / 2 — J7/2 SO splittings in 40 Ca, 56 Ni, and 48 Ca requires (i) a signif- 
icant reduction of the standard isoscalar spin-orbit strength and (m) strong attractive 
tensor coupling constants. The aim of this paper is to address the consequences of these 
strong attractive tensor and weak spin-orbit fields on total binding energies, two-neutron 
separation energies and nuclear deformability. 



1. Introduction 

Density functional theory (DFT) is a method of choice in large-scale calculations of 
nuclear properties. In spite of certain difficulties related to rigorous formulation of 
the DFT for self-bound systems like atomic nuclei, the method is po tenti ally exact, 



which is guarante ed by the Hohenberg-Kohn-Scham (HKS) theoremjU2J see recent 
discussion in RefsP^ES 

Due to the complexity of nuclear many-body problem in general, and of the 
effective nucleon-nucleon interaction in particular, there exist no universal rules 
for constructing nuclear energy density functional (EDF). In this respect, our in- 
tuition is almost solely based on symmetry properties and practical knowledge, 
accumulated over the years for density-dependent effective interactions of Skyrme^ 
or Gognj0 type being applied within the mean-field (MF) approximation. Free pa- 
rameters of these interactions or functionals are fitted to empirical data. Hence, 
the quality and performance of these methods strongly depends on adopted fitting 
strategies and datasetd^l 

Conventional fitting methods use datasets that are dominated by bulk nuclear 
matter data and by nuclear binding energies of selected double-magic nuclei, with 
essentially no data on the single-particle (s.p.) energies. Such strategies have quite 
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dramatic consequences concerning mostly spin-orbit (SO) and tensor parts of the 
EDF. In particular, they lead to an artificial isoscalar-effective-mass scaling of the 
SO strengths^, contradicting scaling in selected s.p. splittings^, and perpetual 
problems in reproducing evolution of the proton (neutron) s.p. energy splittings 
versus the neutron (proton) shell filling along the isotopic or isotonic chains of 
nuclei. 

The most p romine nt example s of such ch ains in clude: neutron-rich oxyge rP, 
neorP^, sod iurrP ^l^H magnesium^], titaniurr l^^ , or, in medium-mass region, 

I1QI 

antimony^ 1 isotopes. In fact, a non-conventional shell evolution found in these 
neutron-rich nuclei directly motivated the shell-model theorists to introduce the 
so-called monopole shifts, to account for empirical trends. The physical origin of 
these shifts was, in turn, attributed to the shell-model tensor interaction^SZSEI] 
Connection between the monopole-shifts and the tensor interaction was later on 
confirmed within the self-consistent MF ' models ; using ei ther finite-range Gogny 
forced or contact Skyrme interactio d 2 ^ 24 ! 25 ^ ! 27 ! 28 ! 29 ! augmented by a strong 
tensor interaction. 

Single-particle spectra provide for a clear evidence of strong tensor interaction 
and call for a new strategy of fitting the nuclear EDF in gen eral, and of the SO 

IA6IOAI 

and tensor terms in particular, directly to the s.p. dat a 1 I The use of the s.p. 
levels was usually conteste d, becaus e of the isoscalar-effective-mass scaling (m*) of 
s.p. levels. Several authorJ^IElESl argued that the physical density of s.p. levels 
around the Fermi energy can be reinstated only after the inclusion of particle- 
vibration coupling, that is, by going beyond MF. In our opinion, effective EDF 
theories should warrant a proper value of the effective mass through the fit to 
empirical data and readjust other coupling constants to this particular value of m*, 
leading to fairly m* independent predictions. Hence, fitting strategies can include 
information on s.p. levels provided that the s.p. levels are understood through the 
binding-energy differences between doubly -magic cores and the lowest s.p. states in 
odd- A single-particle/hole neig hborJSEEI. Spherical s.p. energies or, more precisely, 
the Kohn-Sham s.p. energies computed in even-even double-magic core should serve 
only as auxiliary quan titi es. 

In our recent studjJ2H, we have proposed a novel fitting strategy of the SO and 
tensor terms in the nuclear EDF. It is based on a direct fit to the ff/2 — /s/2 
SO splittings in spin-saturated isoscalar nucleus Ca, spin-unsaturated isoscalar 
nucleus 56 Ni, and spin-unsaturated isovector nucleus 48 Ca. The procedure allows 
for fixing three out of four coupling constants in this sector, namely, the isoscalar 
strengths of the SO and tensor interactions and the ratio of the isovector coupling 
constants. The procedure indicates a clear need for a major reduction of the SO 
strength and for strong attractive tensor fields. The aim of the present work is 
to address further consequences of strong attractive tensor and weak SO fields on 
binding energies, two-neutron separation energies, and nuclear deformability. 
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Fig. 1. The tt/i 1 1 /2 — ^9t/2 splitting in ^Sn isotopes versus A. Black triangles label empirical data 
taken from Ref. UM, Open and gray symbols represent the SHF results obtained by using the SLy4 
(left) SkO (middle), and SLy5 (right) parameterizations, respectively. Different symbols labeling 
theoretical results follow the SHF minima corresponding to configurations differing in the 
occupancies as indicated in the legend. 



2. Fitting the tensor strengths to single-particle energies 

The s.p. levels constitute one of the main building blocks of the MF method. 
In spite of that, the Skyrme HF (SHF) method that uses forces fitted t o bulk 
nuclear properties performs rather poorly with regard to the s.p. SO splittingJ^HEQ]. 
This is visualized in Fig. [T] showing the 7rl/iii/ 2 — 7iT<77/2 splittings, AE g h, in 
antimony, calculated by using the SLyJSI, SkcP, and SLyP2 parameterizations . 
The SLy4 force strongly overestimates the absolute value of AE g h and fails to 
reproduce the slope of the AE g h(A) curve. The non-standard isovector SO in the 
SkO force helps by reducing, on average, the splitting to the empirical level, but 
does not change the slope of the AE g h(A) curve. Finally, in SLy5, the inclusion 
of tensor terms changes the slope, but shifts the theoretical curves in a wrong 
direction. The latter observation suggests that the fit to masses leads to values of 
tensor coupling consta nts tha t are at variance with those deduced from the s.p. 
level analysis, see RefsP^H. However, one should point out that the nlhn/ 2 — 
irlgr/2 splittings depend upon many factors including, apart from the SO and tensor 
fields, the effective mass, centroid energies of the £ — 4 and I — 5 sub-shells, 
and strong polarization effects. Hence, conclusions concerning the SO and tensor 
coupling constants that are deduced solely from these data should be considered to 
be tentative. 

It is well known, see Refs.HIE^l, that the tensor interaction strongly modifies 
the SO one-body potential. In the spherical-symmetry limit, the isoscalar (t = 0) 
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Fig. 2. The neutron (top) and proton (bottom) I/7/2 — I/5/2 SO splittings in 40 Ca, 48 Ca, and 
56 Ni. Black symbols show the mean empirical values taken from 

Refsl38 39 

and open dots denote 

the SkO results. Open diamonds represent the results obtained by using the SkO jv functional of 
which includes strong attractive tensor terms and a reduced SO strength. 



and isovector (t = 1) SO one-body potentials read: 



1 

27 



Cl'Jt(r) - C, 



yjdpt 
dr 



L S, 



(1) 



where C/ and C t VJ are the tensorial and spin-orbit coupling constants, see for 
example RefP. The tensor field depends upon the radial component of the spin- 
orbit vector density J t = - Jt(r) that measures the spin-asymmetry of the nucleus 
and can rapidly vary with particle numbers. On the contrary, the second term in 
Eq. ([T]) , which is due to the conventional two-body spin-orbit interaction, depends 
on the radial derivative of the particle density pt, which varies relatively slowly 
with particle numbers. Such a contrasting behavior of the two major constitu ents 
of the SO potential can be actually used to fit the coupling constants to dataP^l. 
The idea is visualized in Fig. [21 which shows the 1/7/2 — I/5/2 SO splittings in 
40 Ca, 48 Ca, and 56 Ni. These splittings form a very distinct pattern that cannot be 
reproduced based solely on the conventional SO potential. Indeed, the I/7/2 — 1/5/2 
SO splittings in 40 Ca, 48 Ca, and 56 Ni are fairly constant when calculated using, 
for example the SkO force, see curve marked by open dots in Fig. [2] It reflects the 
fact that the neutron and proton radial form-factors ^ almost do not change when 
going from 40 Ca through 48 Ca to 56 Ni. At the same time the neutron and proton SO 
vector densities J(r) change rapidly when going from the isoscalar spin-saturated 
40 Ca to the isoscalar spin-unsaturated nucleus 56 Ni, and, finally, to the isovector 
spin- unsaturated nucleus Ca. This allows for a simple and intuitive three-step 
fitting procedure^! f the C^ J in 40 Ca, C J in 56 Ni, and C(/CY J ratio in 48 Ca. 
This procedure leads to (i) a significant reduction in the isoscalar SO strength and 

(ii) strong attractive tensor coupling constants. It systematically improves such s.p. 

1901 

properties as the SO splittings and magic-gap energies^ 21 , but leads to deteriorated 
nuclear binding energies. 



3. Tensor interaction and the nuclear binding energies 

As discussed in RefP^l, the tensor contribution to the nuclear binding energy 
shows interesting generic topological patterns closely resembling those of the shell- 
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Fig. 3. Tensor contributions to the total binding energy c alcu lated by using the spherical Hartree- 
Fock-Bogoliubov model with the SLv4t fun ctional of Ref|2§l Vertical and horizontal lines indicate 
the tensorial magic numbers. From RefJSl 
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Fig. 4. Differences between theoretical and experimental binding energies (positive) in selected 
double-magic nuclei. Black dots represent the results obtained using conventional SLy4 force. 
Black and white diamonds label the results calculated using the SLy4y and the SkOjv functionals, 
respectively. 



correction, see Fig. [3] The single-particle tensorial magic numbers at N(Z) = 14, 32, 
56, or 90, corresponding to the maximum spin-asymmetry in the l/7/2©2p 3 / 2 , 
l59/2®2ri 5 /2 and l/iii/2©2/ 7 / 2 spherical s.p. configurations, respectively, are clearly 
seen in the figure. Note, that the calculated tensorial magic numbers are shifted 
due to configuration mixing toward the classic magic numbers of N(Z)=8, 20, 28, 
50, and 82. The topological features shown in Fig. [3] are fairly independent of a 
specific parameterization of the force. Indeed, they simply reflect the order of s.p. 
levels, which is rather unambiguously established and relatively well reproduced by 
the state-of-the art nuclear MF models, at least in light and medium-mass nuclei. 

Values of the tensor and SO strengths de duced from the s.p. properties are 
at variance with those obtained from mass A large reduction of the SO 



November 7, 2008 19:32 WSPC/INSTRUCTION FILE kazi2008-07 



6 W. Satuia, J. Dobaczewski, P. Olbratowski, M. Rafalski, T.R. Werner, R.A. Wyss, and M. Zalewski 



18 20 22 24 26 28 30 



> 10 

o 

S 5 



> 1 

01 

§ 



CO 



oxygen 




^TAME03 



18 20 22 24 26 28 30 

A 



Fig. 5. Two-neutron separation energy S2n (upper part) in oxygen isotopes. Empirical data S3] 
are labeled by black dots. Theoretical values obtained using the SkO and SkOy/ functionals are 
marked by gray and white dots, respectively. Lower part shows contribution of the tensor term to 
the S2n as a function of shell filling. 



strength, which is particularly strong for the low-m* forces like SLy4, has a partic- 
ularly destructive impact on theoretical binding energies, see Fig. [4] 

A multidimensional fit to masses shows that the mass performance of the SLy4r 
force can be improved by a tiny refinements of the remaining coupling constants^ 21 , 
however, the accuracy of the original SLy4 cannot be regained. This indicates that 
a spectroscopic quality parameterization that would perform reasonably well on 
binding energies must have lar ge e ffective mass, m* > 0.9. One of the candidates 
is the SkOy' functional of RefM This functional, at least for the classical set 
of double-magic nuclei shown in Fig. 31 is of a similar accuracy as SLy4, and it 
outperforms both the SLy4y and SLy4<r m m of RefP. 

This allows for reasonable quantitative estimates of the tensor influence, for 
example, on two-neutron separation energies, potential energy surfaces (PES's), 
and onset of deformation. An example of calculation of the two-neutron separation 
energies for oxygen nuclei is shown in Fig. [5l One clearly sees here the way the 
tensor interaction induces in oxygen isotopes a breaking of stability against the 
two-neutron emission around 26 0. Indeed, as shown in the lower panel of the figure, 
in 26 O a decrease of S%n is directly related to the d 3 / 2 sub-shell occupation that 
reduces the spin- asymmetry and tensor contribution to the binding energy. 

By deforming the nucleus one can easily change the spin asymmetry and, in turn, 
tensor effects. Fig. [5] shows the PES's versus quadrupolc deformation in 80 Zr (left) 
and 120 Sn (right), calculated by using the quadrupole-constrained HFB method. 
At the spherical shape, nucleus 80 Zr is spin-saturated. By deforming the system, 
one increases the spin-asymmetry by enforcing the occupation of the g g / 2 sub- 
shell. By adding to SkO a strong attractive tensor field (SkOrx), one pulls the 
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Fig. 6. Potential energies versus axial quadrupole deformation parameter /?2 in 80 Zr (left) and 
120 Sn (right), calculated by using the SkO, SkOyxi SkO-jv, and SLy4 functionals. See text for 
details. 

deformed minimum down. The consecutive reduction of the SO strength (SkOT') 
provides for a compensation mechanism, and it shifts the g g / 2 sub-shell and the 
deformed minimum up in energy, back to its original position obtained for the 
SkO functional. In 56 Ni, a similar compensation mechanism is found in the yrast 
super-deformed bands, see 

RefP. In case of 120 Sn, the PES's calculated by using 
the SkO and SkOr' parameterizations are again very close to each other. Note 
however, that both these curves differ substantially from the PES calculated using 
the SLy4 parameterization. 



4. Summary 

In this study, we discussed specific nuclear-structure effects induced by strong at- 
tractive tensor fields and weak spin-orbit field, which result from direct fits of the 
coupling constants to the SO splittings. In particular, we showed that contributions 
to the nuclear binding energies that are due to the tensor field show a generic magic 
structure with tensorial magic numbers at iV(Z)=14, 32, 56, or 90 corresponding to 
maximum spin-asymmetry in I/7/2 © l<?9/2 © 2d 5 / 2 and I/2.11/2 © 2/7/2 

single-particle configurations. We also demonstrated that it is possible to construct 
a functional being able to simultaneously reproduce the /s/2 — /V/2 SO splittings in 
40 Ca, 56 Ni, and 48 Ca nuclei and binding energies of doubly magic spherical nuclei. 
By using this particular functional, we discussed specific structural effects, pertain- 
ing to strong tensor terms, exerted on two-neutron separation energies in oxygen 
isotopes and PES's in 80 Zr and 120 Sn. In particular, in the context of nuclear defor- 
mation properties, we discussed compensation mechanism between the attractive 
tensor fields and weak SO field. 

This work was supported in part by the Polish Ministry of Science under Con- 
tract No. N N202 328234, by the Academy of Finland and University of Jyvaskyla, 
within the FIDIPRO programme, and by the Swedish Research Council. 
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